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! Abstract 

^ ' The nuclear effects in the spm-dependent structure function gi of the deuteron 

. are studied in the kinematics of future experiments at CEBAF, {u < 3 GeV, < 

I 2 GeV^). The magnitude of the nuclear effects is found to be significantly larger 

^ ' than the one occurring in deep inelastic scattering —f oo, — > oo). We discuss 

the mechanism leading to large effects in the region of the nucleon resonances. A 
possibility to measure the neutron structure functions in the CEBAF experiments 
with deuterium is analysed, and conclusions about the experimental study of the 



^ ' dependence of the Gerasimov-Drell-Hearn Sum Rule for the neutron are drawn. 

I. Recently it has been proposed [Q at CEBAF to experimentally study the spin- 
^ ' dependent structure function (SF) of the neutron, g^, in a wide interval of energy, u (0.2 — 

■ 3 GeV), and momentum transfers, (0.15 — 2 GeV"^), using polarized deuterium and 
^He targets. These experiments will shed light on a number of quantum chromodynamics 
(QCD) sum rules and will help establish a connection between results predicted by low 
energy theorems (Q^ — *■ 0) and perturbative QCD (Q^ ^ m^, m being the nucleon mass). 

Of particular interest is the dependence of the Gerasimov-Drell-Hearn sum rules |^ 
for the neutron and proton, and the connection of those sum rules with the Bjorken and 
Ellis- Jaffe sum rules (see review ^ and references therein). 

Keeping in mind the lessons of the EMC-effect, one might expect that nuclear cor- 
rections could play an important role in estimating the neutron SF from the combined 
nuclear and proton data |^ ^. In the region of finite ~ and z/ ~ m, nuclear correc- 
tions are much more important than those in the deep inelastic limit [^0- In this talk 
the role of nuclear structure effects in electron-deuteron scattering in the resonance region 



"'^Talk presented at the 9'^ Amsterdam miniconference "Electromagnetic studies of the deuteron", 
NIKHEF, Amsterdam, February 1-2, 1996. 
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will be discussed, paying special attention to the procedure of extracting the neutron SF 
from the deuteron data in the kinematics of future experiments at CEBAF. 

II. The nucleon contribution to the deuteron structure functions is usually calculated 
by weighting the amplitude of electron scattering on the nucleon with the wave function of 
the nucleon in the deuteron (for recent developments see e.g. 0, |Tl| and references 
therein). For the spin-dependent SF the most important effects are those of the Fermi 
motion and depolarizing effect of the D-wave. Additional effects, such as off-mass-shell 
effects or nucleon deformation, are found to be small |jT2|, |13[. For finite values of and 
I/, the deuteron SF (jff (x, Q^) reads as follows: 



(2) 



ymax{X')Q ) 



Vmin i^^^Q^ ) 



y 



where gi = {gi+g^)/2 is the isoscalar nucleon SF and \E'^^(fc) the deuteron wave function 
with spin projection M. In the rest-frame of the deuteron, with q opposite the z-axis, 
kinematical variables are defined as: 



kq = iy{kQ + ^(x, Q'^)k'j), k^ = m + to - /2m 



X 



QV^kq, 



(3) 
(4) 



where = —2.2246 MeV is the deuteron binding energy. The limits ymin{max){x,Q^) 
are defined to provide an integration over the physical region of momentum in and 
to take into account the pion production threshold in the virtual photon-virtual nucleon 
scattering^. Since both ymin{x, Q"^) and ymax{x, Q"^) are solutions of a transcendent equa- 
tion, explicit expressions for them cannot be given. However, in our numeric calculations 
they are accurately taken into account. 



Eqs. (|T])-(|D have the correct limit in the deep inelastic kinematics [Q^ — > oo. 



oo). In this case: ^(x, Q^) 
formula for the deuteron SF 



yn 



~^ yrr 

nH is recovered: 



M£)/m^ and the usual convolution 



Mo/m 



g^{x,Q' 



-9n^/y,Q')fDiy). 



(5) 



Equation (^) defines the spin-dependent "effective distribution of the nucleons", fo, which 
describes the bulk of the nuclear effects in g^. The main features of the distribution 
function, foiy), are a sharp maximum aX y = 1 + to /2m 0.999 and a normalization 
given by (1 — 3/2Pd) (Pq being the weight of the D-wave in the deuteron). As a result, in 
the region of medium values of x ~ 0.2 — 0.6, the deuteron SF (^f (x) is slightly suppressed 



^For X not too close to the limit of single-nucleon kinematics, x 
can be disregarded 



1, the quasi elastic contribution 
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by a depolarization factor, (1— 3/2Pd) xg^[x), compared to the free nucleon SF. However, 
the magnitude of this suppression is small (~ 1%) and this is why it is phenomenologically 
acceptable to extract the neutron SF from the deuteron and proton data by making use 
of the following approximate formula: 

g^ix, Q') ^ (l - ^Pz.) (g^x, Q') + 9{{x. Q'))/2. (6) 

In addition, when integrated over x, eqs. (||) and (||) give exactly the same results (F = 
Jdxgi{x)), i.e. 

= (l- Ipd) (r„(Q2) + F,(g2))/2, (7) 
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which allows one to define exactly the integral of the neutron SF F„ from the deuteron 
and proton integrals, without solving (||). 

Eqs. (|l|)-(||) at finite values of and u are more sophisticated than the corresponding 
equations in the deep inelastic limit. In particular, they do not represent a "convolution 
formula" in the usual sense, since the effective distribution function fo and the integration 
limits are also functions of x. This circumstance immediately leads to the conclusion that, 
in principle, when integrals of the SF are considered, the effective distribution can not be 
integrated out to get the factor similar to (1 — 3/2Pd) in (^). Another interesting feature 
of formulae (|l])-(|D is the Q^-dependence of fn and ymin,(max){x,Q'^). If we again limit 
ourselves to the discussion of the integrals of SF, one concludes that the Q^-dependence 
of such an integral is governed by both the QCD-evolution of the nucleon SF and the 
kinematical Q^-dependence of the effective distribution of nucleons. 

Thus, we have established that in the non-asymptotic regime, equation (J^), in princi- 
ple, does not hold. Furthermore, it is not clear whether an equation similar to could be 
applied in this region. Indeed, we are discussing the kinematical conditions pertaining to 
nucleon resonances, where the "elementary" nucleon SF explicitly exhibits Breit-Wigner 
resonance structures corresponding to the excitations of the nucleon by the photon and 
one expects that the Fermi motion and binding of nucleons will result in a shift and 
smearing of the resonance structures. However, one can hope the actual effects will be 
quantitatively small so that eqs. similar to (|^) and (|^) can phenomenologically still be 
valid. 

III. In our numerical estimates we use a reliable parametrization of the proton and 
neutron SF given by Burkert |]T3, which takes into account several nucleon excitations and 



provides a reasonable description of the available nucleon data in the resonance region. 



Using the Bonn potential model for the deuteron wave function |]T5[, we carry out a 
realistic calculation of the deuteron SF, gi{x, Q"^) in the region of nucleon resonances. 

In Fig. 1 and 2 the results of the calculation of the deuteron SF, g^i^x^Q"^) and 
F^{x,Q'^) at = 0.1 GeV and 1.0 GeV, are compared with the input of the calcu- 
lation, i.e. the isoscalar nucleon SF, g^{x,Q'^) and F2{x,Q'^). It can be seen that the 
role of nuclear effects in the resonance region is much larger (up to ~ 50% in the maxima 
of the resonances), than in the deep inelastic regime (~ 7 — 9%, depending upon the 
models [H, 0, 13, O], resulting in ~ 6 — 7% from the depolarization factor (1 — 3/2Pd) 



3 



and ~ 1 — 2% from the binding effects and Fermi motion). Such a drastic effect is a 
consequence of the presence of the narrow resonance peaks in the nucleon SF. 

Indeed, let us write the elementary nucleon SF as a sum of smooth background con- 
tribution, Qi'''^, and several resonances, (pi, which can be both positive and negative: 

g^x, Q') = 9r''%x, Q') + ^ Q^). (8) 

i 

The size of the effect of the Fermi motion and binding, for the smooth function, gi'^^ , 
is similar to the one in the deep inelastic regime. The smearing of the nucleon resonances 
in the deuteron SF is estimated by formula: 

x/(xi-Ai) 

0(x,g2) = 0,(x„Q2)__ / dyf{y,Xi,Q'), (9) 



<y > 

x/{xi+Ai) 

where Xi and 2Aj are the position and width of the i-th resonance, < y >~ 1. In 
deriving formula (^) we approximated the resonances by the rectangles of height (f){xi) 
and width 2Aj. This estimation shows that resonance is smeared over the entire region 
of X and is strongly suppressed everywhere, if limits of integrations in r.h.s. of eq. (^) 
are close or, the same, if 2Aj is small. Formula (§) predicts that resonances are more 
suppressed if Aj is smaller and if Xi is larger. We present in Fig. 3 the results of pedagogical 
calculations, aiming to illustrate the features of the formula (|^). One can explicitly see 
that narrow resonance structures at high x are strongly suppressed by the convolution 
(see also behaviour of resonances on Figs. 1 and 2). 

Fig. 4 shows the results of the extraction of the neutron SF from the deuteron and 
proton data by using the approximate formula @, which we believe to give an upper 
limit of the possible errors in this extraction. To emphasize the role of nuclear effects in 
the region of finite Q^, the extracted neutron SF is compared with the original (input in 
the calculation) parametrization of the neutron SF. The use of the approximate formula 
(§) appears to be in some regions completely unreliable. This can be easily understood 
as follows: the proton and neutron SF have similar behavior in the resonance region, in 
that the positions of the nucleon resonances are the same for both of them, whereas the 
resonances in the resulting deuteron SF are smeared and shifted, compared to the isoscalar 
SF. Therefore, the subtraction of the proton SF from the deuteron one, in the maximum 
of the former, can result in a minimum for the neutron SF, instead of a maximum. The 
conclusion of our analysis is that nuclear effects in the resonance region are very specific 
and the approximate formula (^ does not work, even for the crude extraction of the 
neutron SF. Obviously, another method of extracting the neutron SF should be used. 



IV. In ref. [0 a rigorous method of solving eq. (|^) for the unknown neutron SF 



has been proposed and apphed in the deep inelastic region. It has been shown that this 
method, which works for both spin-independent and spin-dependent SF's, in principle 
allows one to extract the neutron SF exactly, requiring only the analyticity of the SF. 
It can also be applied by a minor modification to the extraction of the SF at finite Q^, 
which is our present aim. 
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The basic idea is to replace the integral equation by a set of linear algebraic 
equations, KQ^^ = Go, where is a square matrix (depending upon the deuteron model), 
Qd is a vector of the experimentally known deuteron SF and is a vector of an unknown 
solution. Changing the integration variable in (^), r = x/y, we get 



9i{x,Q')= J dTg^{T,Q')^Ux/T,ax,Q')), (10) 



where TYnm^X^Q ) x/t/max^XjQ ), T^dxi^XjQ ) Xfnaxi^Q ^/Umin^X^Q ) and Xmaxi^Q ) is 

defined by the pion production threshold in virtual photon-nucleon scattering. Let us 
assume that the deuteron SF has been measured experimentally in the interval {xi,X2) 
and a reasonable parametrization for the SF is found in this interval. Then, dividing both 
intervals {xi, X2) and (rmm, Tmax) iiito small parts, one may write: 

9?ix,, Q') ^ E 9i{r„ Q') I -fnix^/T, Q')dr, t = l...N, (11) 



where tj = Tmin + h{j — 1/2) and h = {Tmax — Tmin)/N. Equation ([Til) is already explicitly 
of the form Qd = KQn, therefore the usual linear algebra methods can be applied to solve 
it. 

Note that the range of variation of r is larger than the one for x. Therefore, in principle, 
the SF of the deuteron, experimentally known in the interval (xi, X2), contains information 
about neutron SF in wider interval (for example, in deep inelastic regime r^m ~ x/2 and 
^max = !)• However, extracting information beyond the interval fmin = Xi to fmax = X2 
is almost impossible in view of the structure of the kernel of eq. ( pUj ) and the kinematical 
condition of planned experimental data ||16[. We have to redefine the kernel of eq. ([T0|) to 



incorporate new limits of integration fmin = Xi and f^ax = X2 ||T6|. 

The procedure of solving the eq. (H) in the kinematical region of finite and u will 
be presented elsewhere in details; here we only stress that the method works with good 
accuracy. To check it, we calculated the deuteron SF by formula (^ with the nucleon 
SF g^lXjQ"^) from ref. [|l^ and the deuteron wave function of the Bonn potential jTS 



Then the obtained gf{x, Q^) has been used as "experimental" data to calculate the vector 
in (0); the matrix K has been calculated by using the same deuteron wave function. 
Equation (|TTp has been solved numerically for various "experimental" situations (changing 
the "measured" interval (xi,X2), for different Q^, etc.). The obtained solution, i.e. the 
extracted neutron SF, has been compared point by point with the input to the calculation 
of gi- We found that method is stable and allows one to unfold the neutron SF with 
errors not larger than 10~^, which is much smaller than the expected experimental errors. 
(Note, that all results and conclusions are valid for both polarized and unpolarized SF.) 

V. In this section the role of nuclear corrections in the analysis of the integrals of the 
SF, such as the GDH Sum Rule will be discussed. A very important observation has been 
made in the deep inelastic limit, the exact formula (|^) and the approximate formula (^) 
give the same result for the integral of the neutron structure function, g^^x^Q"^ ^ m^) 
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(see eq. (|^). Likewise, we consider a possibility to apply the eq. (||) to the estimate of the 
integral of the nucleon SF from the deuteron data. The applicability of the approximate 
formula in the deep inelastic region is based on the conservation of the norm of the 
distribution, f{y), by the convolution formula. This circumstance can not be immediately 
extended to the case of the resonanse region, (i) the convolution is broken in eq. (|^) and 
(ii) the normalization of the function f{y,x,Q'^) is different from one of f{y). However, 
the size of the effects is not too large and they should not lead to large errors if we use 
eq. (0). 

In order to understand the deviations of the integral of the deuteron SF in the reso- 
nances region from eq. (|^, let us evaluate the deuteron SF, (||), using the presentation 
(|) of the nucleon SF: 

g^{x,Q') ^ g^^^\x/{y),Q') Uy,x,Q')dy + Y.M^,Q'), (12) 

where the first term in the r.h.s. of eq. (0) is obtained using an expansion at the point of 
sharp maximum of distribution function, /, and the second term is defined by (^. Next, 
we define an auxiliary function rig// as: 

Vmax iS^iQ ) 

ne/y(x,Q2) ^ (^1 _ Ipj^^ J f^(y,x,Q')dy. (13) 

Thus, it can be seen that Ueff represents the "effective number" of nucleons "seen" by 
the virtual photon in the process when the virtual photon is absorbed by the nucleon and 
at least one pion is produced in the final state. Obviously, in the deep inelastic limit, 
rieff = 1 corresponds to the normalization of the deuteron wave function. The effect of 

on rieff consists in narrowing the interval on x from to XmaxiQ"^)- At finite Q^, 
Ueff is close to 1, in so far as x (or u) is far from the threshold, and rze// rapidly falls 
down only very close to the threshold. Figure 5 illustrates the behaviour of neff{x,Q'^) 
as a function of x. The arrows indicate the kinematical limits Xmax at a given Q^. The 
function neff{x, Q^) significantly differs from unity only when x — > x^ax where we expect 
the nucleon SF to be rather small (see Figs. 1,2). 

Then, integrating ( [T^ ) we get: 

{Xtr ^ 
J rfx(7f''^(x,g2)ne//(x,g2) +^2A0,(x„Q2)ne//(x„g2)| .(14) 

Noting that the integral of the nucleon SF, (^, is 

r^(Q2) = J dxg^^'''{x,Q') + Y.2A<p,{x„Q'), (15) 

* 

and neff[x) is close to 1, we expect the quantity in curly brackets in the r.h.s. of eq. ([T^) 
is close to the integral of the nucleon SF, (|T5|). 
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Therefore, we expect only small corrections to the integral of the deuteron SF com- 
pared to the (^). This effect can be accounted for by a new equation: 

Td{Q') = (i - \pd) iVe//(g^)(r„(g2) + t,{q'))/2, (i6) 

where Neff{Q^) ^ J dxneff{x,Q'^); eq. (p!6| ) and the integral of eq. (0) represent the 
definition of the effective number Ne/flQ"^), which depends upon the form of the nucleon 
SF, g^; and, since this is expected to strongly oscillates (see Fig. 1), even the sign of the 
correction can vary. For instance, we obtain using the SF from ||14|| , 

NeffiQ^ = 0.1 GeV^) = 1.02, N^jfiQ^ = 1.0 GeV^) = 0.997, (17) 

i.e. a rather small effect (+2% and —0.3% correspondly). Therefore eq. (|T6|) appears to 
be reliable for estimating the integrals of the SF: setting NeffiQ"^) = 1 does not lead to 
errors larger than 3% for = 0.1 — 2.0 GeV^. 

VI. In conclusion, we have shown that the effects of nuclear structure in the extraction 
of the neutron SF in the resonance region are much more important than in the deep 
inelastic scattering. We have explained how the correct neutron SF can be firmly extracted 
from the combined deuteron and proton data. At the same time, we have found that the 
integrals of the SF, such as the GDH Sum Rule, can be estimated with accuracy better 
than 3% by the simple formula (|^) which is also valid in deep inelastic region. 
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Figure 1. The spin dependent structure functions gi{x,Q'^) for two values of 
Q^. The deuteron SF (sohd hne) is compared with the isoscalar nucleon SF 
(dotted hne) used as input into the calculation in eq. (Ef). 
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Figure 2. The spin independent structure functions F2{x,Q'^) for two values of Q^. 
The deuteron SF (sohd hne) is compared with the isoscalar nucleon SF (dotted hne) 
used as input into the calculation in similar to eq. (0). 
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Figure 3. Illustration for formula (|]). Examples of calculation of the convolution 
(BI), solid lines, with narrow resonance structures, dashed lines. 



11 



0.1 
0.0 - 

-0.1 - 
-0.2 - 
-0.3 - 

-0.4 - 
-0.5 



-0.6 



g,(x,q''=0.1 GevVe'^) 




Neutron (extracted) 
Neutron (input) 



0.00 - 



-0.05 



g,(x,Q'=l GevVc") 



-0.10 




0.0 0.1 0.2 

X 



— I -0.15 - 
0.3 0.{ 



Neutronf extracted) 
Neutron(input) 



0.2 0.4 0.6 0.[ 

X 



Figure 4. The extracted neutron SF (dotted line) by approximate formula 
eq. in comparison with the original parametrization (solid line) used into 
the convolution formula (0). 
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neff(x,Q ) 
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Figure 5. The effective number neff{x,Q'^) characterizing the additional nu- 
clear effects in the GDH integral, eq. (|T3|). 
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